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ABSTRACT 

The  combustion  or  afterburning  of  fuel-rich  rocket  exhaust  with  the  atmosphere  may  result  in 
large  infrared  radiation  emissions  which  can  play  a  significant  role  in  the  design  of  missile  base 
components  and  missile  defense  systems.  Current  engineering  level  models  neglect  turbulence- 
chemistry  interactions  and  typically  underpredict  the  intensity  of  plume  afterburning  and  afterburning 
burnout.  To  evaluate  the  impact  of  turbulence-chemistry  interactions,  an  assumed  pdf  model  was 
applied  to  missile  plume  simulations  of  a  generic  booster.  Simulation  results  reveal  turbulence- 
chemistry  interactions  to  have  a  large  impact  on  plume  signatures  as  afterburning  burnout  was 
approached. 
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NOMENCLATURE 

model  calibration  coefficients 

mixture  specific  heat  at  constant  volume 

mean  mixture  specific  heat  neglecting  turbulent  fluctuatioiis 

mixture  specific  internal  energy 
Ath  species  specific  internal  energy 
specific  internal  energy  variance 
mixture  specific  enthalpy 
turbulent  kinetic  energy 
pressure 

single-point  joint  probability  density  function  (pdf)  of  density, 

temperature  and  species  mass  fraction 

marginal  pdf  of  scalar  energy 

marginal  pdf  of  temperature 

marginal  pdf  of  density 

Prandtl  number 

turbulent  Prandtl  number 

scalar  energy 

universal  gas  constant 

sum  of  the  species  mass  fraction  variance 

Schmidt  number 

turbulent  Schmidt  number 

time 
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A 

T 

= 

temperature 

Train,  Tmax 

= 

Pt  pdf  minimum  and  maximum  realizable  temperatures 

T0,  Tj 

Pj  pdf  minimum  and  maximum  absolute  realizable  temperatures 

<  T"1  > 

= 

temperature  variance 

Uj 

= 

yth  velocity  component 

wk 

= 

Mi  species  reaction  rate 

wk 

= 

Mi  species  molecular  weight 

Xj 

- 

yth  spatial  coordinate 

Xrefl 

distance  between  missile  nozzle  exit  plane  and  the  plume  barrel  shock 

reflection  point  at  the  axis  of  symmetry 

Yk 

Mi  species  mass  fraction 

Pi,  @2,  Pk 

= 

beta  function  pdf  parameters 

5 

= 

plume  transverse  length  scale 

A 

= 

delta  function 

s 

" 

turbulent  kinetic  energy  dissipation  rate 

r 

ratio  of  mixture  specific  heats 

r 

= 

gamma  function 

V 

kinematic  viscosity 

vT 

= 

turbulent  eddy  viscosity 

p 

= 

density 

e 

= 

normalized  temperature 

Superscripts 

If 

= 

Favre  fluctuating  component 
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A 


Other  marks 


=  time-averaged  quantity 

<  >  =  Favre- averaged  quantity 

I.  INTRODUCTION 

The  afterburning  of  missile  exhaust  with  the  atmosphere  may  result  in  large  infrared  radiation 
emission  which  can  be  a  major  contributor  to  the  heat  transfer  to  the  missile  base.1  As  altitude  is 
increased,  the  radiative  component  of  the  heat  transfer  rate  will  eventually  show  a.  large  drop  in 
magnitude  resulting  from  the  shutdown  or  cessation  of  afterburning  in  the  plume  (see  Fig.  1). 
Shutdown  in  this  context  does  not  refer  to  the  termination  of  the  missile  engine,  but  to  the  cessation 
of  the  combustion  taking  place  between  the  missile  exhaust  and  the  atmosphere,  occurring  with 
continuous  engine  operation.  The  total  amount  of  time  this  afterburning  occurs  and  its  rate  of  decay 
during  the  shutdown  event  will  determine  the  total  heat  transfer  to  the  body  and  establish  a  design 
criteria  for  the  components  in  the  base  region.  The  radiative  emission  during  the  shutdown  event 
also  has  important  implications  for  the  development  of  missile  defense  systems.  Consequently,  a 
need  exists  to  accurately  characterize  missile  plume  afterburning  shutdown  or  cessation  events. 

The  character  of  afterburning  shutdown  events  has  been  observed  to  vary  among  different 
propulsion  systems.  For  example,  during  shutdown  some  systems  exhibit  a  very  rapid  drop  in  total 
radiant  intensity  over  a  narrow  altitude  range,  whereas  others  shutdown  gradually  over  a  wide  range. 
The  prediction  of  afterburning  shutdown  using  engineering  level  modeling  techniques3  has  been 
relatively  successful  for  systems  that  exhibit  the  gradual  drop-off  type  of  shutdown  event.  However, 
these  techniques  have  been  less  successful  in  predicting  shutdown  for  systems  that  exhibit  the  rapid 
drop-off  b  ehavior . 

To  improve  understanding  into  the  physical  mechanisms  producing  afterburning  shutdown, 
Calhoon4  investigated  the  shutdown  characteristics  of  a  generic  amine  booster  system  within  the 
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framework  of  a  computational  parametric  study.  Several  mechanisms  were  investigated  to  explain 
the  shutdown  behavior  of  this  system  which  included:  1)  a  plume  shear  layer  relaminarization 
phenomenon,  2)  a  Damkohler  number  effect  and  3)  a  classical  flame  extinction  mechanism.  The 
relaminarization  mechanism  was  found  to  be  implausible  because  plume  shear  layer  Reynolds 
numbers  were  well  above  transition  after  shutdown  had  occurred.  The  Damkohler  number 
mechanism  was  found  to  be  plausible  and  indeed  the  only  shutdown  mechanism  modeled  within 
most  commercially  available  codes.  This  is  a  consequence  of  the  assumptions  used  to  model  the 
mean  reaction  rate  in  the  species  conservation  equations.  The  Damkohler  number  mechanism  was 
found  to  be  responsible  for  the  gradual  drop-off  rates  of  plume  radiation  produced  by  most  codes. 
The  classical  flame  extinction  mechanism  was  found  to  be  a  possible  explanation  for  the  observed 
rapid  shutdown  behavior  of  some  systems. 

The  flame  extinction  phenomenon  investigated  by  Calhoon4  is  a  result  of  the  interaction  of 
turbulence  and  chemistry  at  small  spatial  scales.  The  model  used  by  Calhoon4  to  account  for  this 
phenomenon  was  an  elementary  one  applicable  to  high  speed  compressible  flows.  Other  models  do 
exist  to  more  generally  account  for  turbulence-chemistry  interactions  that  are  applicable  to 
compressible  flows.  Among  the  most  comprehensive  techniques  are  the  assumed  probability- 
density-function  (pdf)  method5  and  the  compressible  extension  of  the  pdf  evolution  equation 
method  ’  .  Though  providing  a  more  accurate  description  of  the  higher  order  statistics  of  the 
turbulent  scalar  fields,  the  pdf  evolution  equation  method  has  not  been  shown  to  yield  significantly 
better  results  than  the  assumed  pdf  method  when  applied  to  compressible  flows.  The  pdf  evolution 
equation  method  is  also  computationally  expensive  and  may  become  intractable  when  applied  to 
flows  including  shock  waves.7  The  assumed  pdf  method,  on  the  other  hand,  is  computationally 
inexpensive  and  offers  a  viable  approach  to  account  for  turbulence-chemistry  interactions  in 
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complex,  large  scale  flows  of  practical  interest.  However,  the  assumed  pdf  method  does  not  directly 
account  for  strain  extinction  effects  and  may  not  accurately  capture  this  phenomenon. 

For  high  speed  flows,  turbulence-chemistry  interactions  have  been  shown  to  enhance  burning 
for  flames  far  from  equilibrium.8  Including  turbulence-chemistry  interactions  within  rocket  plume 
simulations  has  the  potential  to  improve  predictive  capability  because  engineering  level  modeling 
techniques  typically  underpredict  afterburning  plume  temperatures  and  emissions  near  the 
afterburning  shutdown  regime.  Including  turbulent  fluctuations  within  the  radiation  calculations  may 
also  significantly  impact  predictions.9  However,  this  aspect  of  the  problem  was  not  considered  in 
this  study.  The  effect  of  turbulent  fluctuations  was  only  considered  for  the  flowfield  predictions. 

The  objective  of  this  study  was  to  assess  the  impact  of  turbulence-chemistry  interactions  on 
the  afterburning  and  afterburning  shutdown  characteristics  of  a  generic  missile  system.  This 
assessment  will  provide  guidance  for  further  investigation  and  for  the  enhancement  of  engineering 
level  models.  This  work  is  a  continuation  of  a  previous  study4  that  investigated  mechanisms 
influencing  the  afterburning  shutdown  characteristics  of  rocket  exhaust  plumes. 

In  the  following  sections  a  brief  review  of  the  assumed  pdf  turbulence-chemistry  interaction 
model  is  first  given.  The  computational  methodology  used  within  the  simulations  is  then  presented, 
followed  by  a  description  of  the  generic  missile  geometry  and  engine  model  used.  Results  are  then 
presented  for  the  prediction  of  afterburning  shutdown  for  this  missile  configuration  using  the 
assumed  pdf  model  followed  by  conclusions  that  may  be  drawn  from  the  study.  These  results 
indicate  that  turbulence-chemistry  interactions  do  play  a  significant  role  in  the  afterburning  and 
afterburning  shutdown  characteristics  of  missile  exhaust  plumes. 

II.  ASSUMED  PDF  METHOD 

The  effect  of  turbulence-chemistry  interactions  within  the  context  of  Reynolds  averaged 
Navier-Stokes  simulations  (RANS)  appears  as  additional  unclosed  terms  in  the  governing  steady 
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state  conservation  equations.  These  unclosed  terms  include:  1)  the  time  averaged  reaction  rate  term 
T.  in  the  species  conservation  equations  and  2)  temperature-species  correlations  in  the  state 
equations.  The  reaction  rate  is  a  highly  nonlinear  function  of  temperature  and  species  concentrations 
and  its  time  average  may  be  expressed  as, 

00  co  1  11 

*,  =  \\[..\\v<t(p,T,Yl)P(p,T,YkW,...dYKdTdp  (1) 

0  0  0  0  0 

where  T,  p,  Yk  and  wk  are  the  temperature,  density  and  Mi  species  mass  fraction  and  reaction  rate  for 
species  1  <k<K.  P(p,T,Yk )  is  the  single-point  joint  probability  density  function  of  p,  T  and  Yk  and 

represents  the  combined  effects  of  turbulent  transport,  both  large  and  small  scale,  and  molecular 
diffusion.  The  time  averaged  state  equation  for  a  mixture  of  perfect  gases  may  be  written  as, 


K  ^  rpnyn 

+R‘pY—  ‘ 

K  wt 


where  the  overbar  and  brackets  (<  >)  represent  conventional  and  Favre  time  averaging,  respectively. 
Also,  p  is  the  pressure,  R°  is  the  universal  gas  constant  and  Wk  is  the  Mi  species  molecular  weight. 
The  fluctuating  components  T"  and  7/  of  T  and  Yk  are  defined  by  T  =<T>  +T"  and 
Yk  =<  Yk  >  +Yk  ,  respectively. 

Temperature-species  correlations  similar  to  what  is  seen  in  Equ.  (2)  also  appear  in  the  time 
averaged  caloric  equation  of  state,  <h>=<e>+p/p,  relating  the  internal  energy  to  enthalpy. 

Expressing  the  species  specific  enthalpies  in  terms  of  standard  polynomial  curve  fits,  the  time 
averaged  enthalpy  becomes  <h>=  f  (pT ’Y", ...,pT"N Yk)  where  N  is  the  number  of  coefficients  in  the 


To  close  the  governing  flow  equations,  an  approximation  for  the  pdf  P  must  be  specified. 


Within  the  context  of  the  assumed  pdf  approach,  the  shape  of  P  is  specified  in  terms  of  functions 
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which  may  be  parameterized  by  their  lower  moments.  The  first  approximation  used  within  this 
approach  is  to  assume  statistical  independence  of  p,  T and  Yk  so  that  P  may  be  expressed  as, 

P(p,T,Yk)  =  Pf{p)Pr(T)Pg(Yk)  (3) 

where  Pt  and  Pq  are  the  marginal  pdfs  of  p,  T  and  Yk,  respectively.  This  assumption  is 
questionable.  However,  other  assumed  pdf  methods  cast  in  terms  of  a  mixture  fraction  and  chemical 
progress  variable  have  been  shown  to  produce  good  results  for  a  range  of  low  speed  turbulent 
reacting  flow  problems  using  this  assumption.  Consequently,  this  assumption  is  carried  over  to  the 
more  general  formulation  described  here.  This  assumption  also  results  in  <  FT/  >  =  0  in  Equ.  (2) 

since  the  variables  are  not  correlated. 

Baurle  and  Girimaji10  investigated  a  modification  to  this  procedure  which  relaxed  the 
statistical  independence  assumption  and  showed  potential  for  improvement.  However,  further 
research  is  required  to  resolve  realizability  issues  that  were  identified.  Therefore,  this  modified 
approach  was  not  considered  further  in  this  study. 

Consistent  with  work  of  Gaffney,  et  al.5,  the  marginal  pdfs  Pp>  Pt  and  Pq  were  specified  as 
follows.  A  delta  function  was  assumed  for  the  marginal  pdf  of  density  so  that  Pp-  A(p  -  p).  For 


the  temperature  pdf  Pt,  a  beta  function  was  used.  The  beta  function  has  been  shown  to  accurately 
capture  scalar  mixing  in  homogeneous  turbulence. 1 1  The  pdf  Pt  was  taken  as, 
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where  T  is  the  gamma  function.  The  variance  <  0  ”9 "  >  was  explicitly  limited  by  th’e  realizability 
constraint, 


<  9"9"  ><<9  >  (\-  <6  >)  (7) 

The  temperature  limits  Tmin  and  Tmax  were  specified  as, 

T«i»  =max[r„<r>-nV<r"  >] 

I - 

Tm  =  min[r,  ,<  F  >  W<  T"2  >] 

where  the  coefficient  n  was  specified  to  have  a  value  of  3  following  Gaffney,  et  a/.8  The  values  To 
and  Tj  are  problem  dependent  absolute  minimum  and  maximum  realizable  temperature  limits, 
respectively.  Specifying  r„„„  and  Tmax  in  this  manner  gives  the  pdf  a  symmetric  distribution  for  small 
temperature  variance  while  allowing  for  nonsymmetric  shapes  when  the  limits  are  capped  by  To 
and/or  Tj.  Allowing  nonsymmetric  pdf  shapes  also  allows  for  the  physically  realistic  situation  of 

temperature  intensities  being  greater  than  1  (i.e.,  V<  r2>/<r>  i). 

For  the  pdf  Pq,  a  multivariate  beta  distribution  was  used  and  has  the  following  form, 


PQ(Yl,Y2,...,YK)  =  r(A+A+-  +  &)  Yfr\Pl~\..Y/K~x A(1  - T Yk)  (9) 

2  K  T(px)T(P2)...T(PK)  12  '  ^  * 


k=\ 


with, 


Pk  =<Yk  > 


f\ -S  ^ 

\  Q  j 


(10) 


where, 


s=E<n>J,e= E<w> 


(ii) 


k= 1 


k= 1 
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The  quantity  Q  is  the  sum  of  the  species  variances  and  is  called  the  turbulent  scalar  energy.  As  with 
Pt,  a  realizability  condition  must  be  explicitly  imposed.  This  condition  was  Q  <  (1  -  S) . 

With  Equs.  (4)  -  (11),  the  joint  pdf  of  p,  T  and  Yk  is  completely  specified  given  <T>, 
<  T"T"  > ,  <  Yk  >  and  Q.  The  mean  reaction  rate  may  then  be  calculated  from  Equ.  (1).  Numerically 


integrating  Equ.  (1)  can  become  very  computationally  expensive  for  a  large  number  of  species. 
Fortunately,  due  to  the  form  of  pdfs,  integrals  involving  species  and  Pq  can  be  evaluated 
analytically.5  The  integrals  involving  T  and  Pt  were  evaluated  numerically  and  stored  in  look-up 
table  form  as  a  function  of  T and  the  variance  of  T(<T”2  >),  as  done  by  Gerlinger  et  al .  2 

The  remaining  issue  with  respect  to  the  evaluation  of  Equ.  (1)  is  the  specification  of  the 
temperature  variance  and  scalar  energy.  These  quantities  were  evaluated  using  modeled  evolution 
equations  ’  derived  from  the  species  and  energy  equations  and  are  given  by, 


Dp<g>  _  3  -_,V  ,  Vr, <?<£>)  ^__Vr 


=  —  py  (-  +  —)- — - 

Dt  3c j  I  Pr  Prr  dx- 


+  2  Cxpy 


3  <e> 


Prr  I  3c j 


_  s  _  _  3  <u,-> 

■  2C  P  <  g  >  7  -  2(r  -  1)P  <  g  >  — — — 
s  k  3cj 


£pQ = ±\pP- + Zl.)  «  +  2c 

Dt  3c j  I  Sc  ScT  3c j  j  ScT 


S<Yt> 


-2CepQj  +  2j]wtYC 


where  e  is  the  internal  energy,  <g>  is  the  internal  energy  variance,  Uj  is  the  velocity,  y  is  the  ratio  of 
specific  heats  and  k  and  s  are  the  turbulent  kinetic  energy  and  energy  dissipation  rate.  The 
coefficients  v,  v-r,  Sc  and  Pr  are  the  kinematic  viscosity,  turbulent  eddy  viscosity,  and  Schmidt  and 
Prandtl  numbers,  respectively.  The  length  scales  of  the  temperature  and  species  fluctuations  were 
assumed  to  be  proportional  to  the  velocity  length  scale  and  specified  in  terms  of  the  turbulent 
Schmidt  and  Prandtl  numbers,  Scj  and  PrT ,  respectively.  The  last  term  in  Equ.  (12)  is  a  modeled 
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dilatation  term  resulting  from  compressibility  effects.8  The  last  term  in  Equ.  (13)  is  a  chemical 
source  term  which  may  be  evaluated  using  the  assumed  pdfs  as  described  by  Gaffney  et  al.5 

The  temperature  and  the  temperature  variance  were  obtained  from  the  mean  internal  energy, 
<e>,  and  its  variance,  <g>,  following  a  procedure  used  by  Baurle  and  Girimaji..10  First,  the 
temperature  and  species  fluctuations  were  neglected  in  the  definition  of  the  mean  internal  energy  so 

that  the  temperature  may  be  evaluated  in  the  usual  manner  by  solving  <e> 

Second,  to  obtain  <  T"1  > ,  the  specific  heat,  cv,  was  linearized  about  the  mean  temperature  so  that, 

(14) 

G 

where  cv  =cv(<T  >,<  Yk  >). 

The  dissipation  coefficients  Cg  and  Cq  in  Equs.  (12)  and  (13),  respectively,  were  specified  as 
Cg  =  Cq  =  0.544  by  matching  the  modeled  scalar  decay  rate  to  that  obtained  from  the  DNS  results  of 
Eswaran  and  Pope14  and  McMurtry,  et  al 15  for  forced  isotropic  turbulence.  The  production 
coefficients  in  these  equations  were  specified  as  Cj  =0.2  and  C2  =  0.5  as  will  be  discussed  in  a  later 
section. 

III.  COMPUTATIONAL  METHODOLOGY 

The  missile  plume  flowfields  in  this  study  were  analyzed  using  CRAFT16  which  is  a 
structured,  finite- volume  code  that  solves  the  compressible  Navier-Stokes  (NS)  equations.  The 
solver  is  fully  implicit  and  uses  Roe/TVD  numerics  for  the  inviscid  fluxes  and  second-order  central 
differencing  for  viscous  and  diffusive  terms.  CRAFT  is  sufficiently  general  to  simulate  finite-rate 
chemistry  and  multi-phase  flows  and  includes  standard  polynomial  curve  fits  for  the  thermodynamic 
properties.  The  code  includes  a  variety  of  advanced  turbulence  model  formulations.17  However,  in 
this  study,  only  the  standard  compressibility  corrected  k-s  model  was  used.  CRAFT  also  includes  an 
implementation  of  the  assumed  pdf  turbulence-chemistry  interaction  model  described  in  the  previous 
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section.  The  CRAFT  code  achieves  good  computational  efficiency  through  a  fully  parallelized 
implementation  of  the  flow  solver  using  a  combined  shared  and  distributed  memory  (using  MPI)  data 
parallel  model.  Computational  efficiency  is  also  enhanced  using  a  parabolized  form  of  the  NS  flow 
solver  (PNS)  for  parabolic/hyperbolic  regions  of  the  flow.  This  PNS  solver  includes  a  sub-stepping 
option  to  allow  for  axial  grid  refinement  in  regions  where  flow  properties  change  rapidly. 

Radiation  calculations  were  decoupled  from  the  flowfield  simulations  and  were  carried  out  as 
described  by  Ludwig  et  al.n  and  Nelson19.  Within  this  approach,  radiation  transport  equations  were 
solved  using  a  band  model  for  the  gas  phase  absorption/emission.  The  calculations  for  the  present 
study  used  a  wide  band  pass  to  encompass  emissions  from  OH,  CO,  CO2  and  H20.  Total  radiation 
intensity  predictions  were  also  made  using  a  field  of  view  large  enough  to  encompass  the  entire 
plume. 

IV.  ROCKET  PLUME  SIMULATIONS 

The  simulations  considered  in  this  study  were  carried  out  for  a  generic  axisymmetric  amine 
booster  described  in  detail  by  Calhoon4.  The  exhaust  nozzle  was  assumed  to  have  an  area  ratio  of  10 
and  protrude  aft  of  the  booster  base.  The  entire  missile  body-base-plume  was  simulated  using 
CRAFT’s  NS  and  PNS  solvers.  The  elliptic  regions  along  the  body  and  in  the  base  region  were 
computed  using  the  NS  solver.  The  hyperbolic  plume  flow  was  calculated  using  the  PNS  option. 
Each  axisymmetric  solution  used  49,000  grid  points  to  resolve  the  body-base  region  and  700,000 
points  for  the  plume  region.  Each  grid,  both  body-base  and  plume,  was  manually  adapted  around 
regions  with  high  gradients  to  ensure  proper  resolution  of  the  flow  features. 

The  inflow  boundary  conditions  at  the  engine  nozzle  exit  plane  were  specified  from  a 
separate  nozzle  flow  calculation  described  in  detail  by  Calhoon.4  A  standard  9  species  -  10  reaction 
step  mechanism4  for  H2/CO  oxidation  was  used  for  this  calculation  as  well  as  for  the  plume 
simulations.  For  this  nozzle  simulation,  the  flow  was  fairly  uniform  across  the  exit,  except  within 
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the  boundary  layer  region  along  the  wall.  The  nozzle  core  temperature  was  ~  1330  K.  The  fuel-rich 
conditions  of  the  engine  were  evident  from  an  excess  of  CO  and  some  H2  in  the  core  with  O2 
virtually  depleted.  This  nozzle  solution  was  assumed  to  be  invariant  with  respect  to  altitude  and  was 
used  for  all  the  body-base-plume  calculations.  This  solution,  however,  did  not  account  for 
turbulence-chemistry  interactions  within  the  nozzle.  At  this  inflow,  the  temperature  intensity  was 
specified  as  5%  while  the  turbulent  scalar  energy  was  specified  as  0%. 

V.  RESULTS  AND  DISCUSSION 

The  body-base  and  plume  flowfields  were  calculated  at  three  altitudes  (25  km,  30  km  and  35 
km)  both  with  and  without  the  pdf  model  described  earlier.  These  altitudes  were  selected  because 
they  span  the  afterburning  shutdown  regime  for  this  generic  booster,  given  an  assumed  trajectory 
profile.4  The  Reynolds  and  Mach  numbers  based  on  freestream  conditions  and  the  body  radius  were 
Re co  =  4x1 06,  2xl06  and  1x1 06  and  M„  -  2.6,  3.2  and  3.9  for  the  altitudes  of  25,  30  and  35  km, 
respectively. 

Plume  simulations  were  carried  out  at  these  three  altitudes  for  three  different  levels  of 
modeling.  The  first  assumed  turbulence-chemistry  interactions  to  be  negligible.  This  approach  was 
termed  the  ‘laminar’  reaction  rate  model.  The  second  model  included  turbulence-chemistry 
interactions  via  the  temperature  pdf  only,  while  the  third  approach  included  both  the  temperature  and 
species  pdfs.  These  three  levels  of  modeling  were  considered  not  only  to  assess  the  effect  of 
turbulence-chemistry  interactions  on  missile  plume  afterburning  shutdown,  but  also  to  assess  the 
effectiveness  of  the  different  aspects  of  the  assumed  pdf  model  itself.  For  the  pdf  cases,  the  absolute 
temperature  limits,  To  and  Tj,  required  to  evaluate  the  pdf  temperature  limits,  Tmin  and  Tmax,  in  Equ. 
(8)  were  specified  as  follows.  With  the  freestream  temperatures  for  the  specified  flight  conditions 
ranging  from  222  to  239  K,  the  absolute  minimum  temperature  was  specified  as  To  —  210  K.  The 
maximum  absolute  temperature,  Tj,  was  specified  as  2800  K  which  corresponded  to  the  maximum 
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temperature  in  the  boundary  layer  along  the  nozzle  wall.4  With  these  absolute  temperatures,  the 

maximum  realizable  temperature  intensity  was  v<  T"2  >/  <T>~  1.6,  which  corresponds  to  the 
condition  when  freestream  and  nozzle  boundary  layer  fluid  co-exist  in  an  unmixed  state. 

As  described  by  Calhoon,4  the  altitude  conditions  of  25,  30  and  35  km  characterize  the 
afterburning  shutdown  regime  of  the  present  generic  booster  system  when  simulated  using  the 
‘laminar’  reaction  rate  model.  For  this  model,  the  plume  at  25  km  is  burning  vigorously,  while  at  30 
km  it  has  progressed  deep  into  the  shutdown  regime.  At  35  km  the  plume  flame  is  almost  completely 
extinguished.  Application  of  the  assumed  pdf  model  to  these  flight  conditions  was  found  to  have  a 
substantial  impact  on  the  afterburning  characteristics  of  this  system.  For  example,  Fig.  2  presents 
contour  plots  of  temperature  and  CO 2  mole  fraction  at  30  km  for  both  the  laminar  rate  and 
temperature  pdf  models.  Simulations  using  both  the  temperature  and  species  pdf  will  be  discussed 
later.  From  this  figure,  the  plume  is  seen  to  rapidly  expand  due  to  the  highly  underexpanded 
condition  at  the  nozzle  exit.  The  plume  barrel  shock  is  clearly  evident  in  the  temperature  field  along 
with  its  subsequent  reflection  off  of  the  axis  of  symmetry  and  its  interaction  with  the  plume  shear 
layer.  Further  downstream  a  series  of  weaker  reflections  persist  until  the  plume  shear  layer  merges 
with  the  axis.  For  the  laminar  rate  case,  plume  ignition  is  delayed  downstream  until  approximately 
three  times  the  distance  between  the  exit  plane  and  the  barrel  shock  reflection  point.  This  is  evident 
from  the  figure  by  the  delay  in  the  rise  of  CO2  in  the  plume  shear  layer.  This  ignition  delay  also 
results  in  a  long  delay  in  the  plume  shear  layer  temperature  rise.  For  the  temperature  pdf  case, 
however,  it  is  clear  the  model  has  a  significant  impact  on  the  combustion  processes  occurring  within 
the  plume.  For  the  pdf  case,  ignition  occurs  close  to  the  missile  base  as  evidenced  by  the  rapid  rise 
in  shear  layer  C02.  Further  downstream,  the  plume  vigorously  afterbums  resulting  in  higher  plume 
temperatures. 
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The  enhanced  burning  realized  for  the  temperature  pdf  model  at  this  altitude  was  found  for 
the  other  altitudes  as  well.  Fig.  3  presents  a  comparison  of  the  laminar  and  temperature  pdf  rate 
model  predicted  mean  temperatures  along  the  plume  axis  for  the  three  altitudes.  In  this  figure  the 
axial  coordinate  has  been  scaled  by  xreji.  This  length  scale  is  the  distance  between  the  nozzle  exit 
plane  and  the  barrel  shock  reflection  point  at  the  axis  of  symmetry.  This  reference  length  was  found 
,  to  approximately  scale  the  results  with  respect  to  altitude  so  that  a  fixed  value  of  x/xreji  corresponds 
to  the  same  relative  location  in  the  plume  for  each  case.4  From  Fig.  3,  the  pdf  model  is  seen  to 
enhance  afterburning  at  all  altitudes  as  evidenced  from  the  higher  peak  temperatures  in  the  plume 
farfield.  At  25  km,  when  the  laminar  model  is  vigorously  afterburning,  the  difference  between  the 
two  models  is  relatively  small.  However,  as  altitude  is  increased  and  afterburning  shuts  down  for  the 
laminar  model,  the  differences  become  more  pronounced.  At  35  km,  the  pdf  case  is  still  afterburning 
while  the  laminar  rate  case  is  becoming  extinguished  resulting  in  a  temperature  difference  of 
approximately  100  K  in  the  farfield  plume.  As  altitude  is  increased,  the  pdf  model  shows  a  gradual 
migration  of  the  ignition  point  further  downstream  consistent  with  the  laminar  rate  model.  However, 
the  temperature  pdf  model  has  delayed  the  onset  of  afterburning  shutdown  significantly. 

Gaffney,  et  al.&,  found  temperature  fluctuations  to  enhance  chemical  reactions  and  ignition 
for  high  speed  Hj-air  shear  layer  flames  that  were  far  from  equilibrium.  The  same  observation  may 

s. 

be  made  here  with  regard  to  CO /H 2-air  combustion  occurring  in  these  plume  flows.  For  the  present 
plume  application,  the  laminar  rate  model  produces  plume  flames  that  are  far  from  equilibrium  and 
on  the  edge  of  burnout.  When  the  pdf  model  is  applied,  temperature  fluctuations  enhance  product 
formation  resulting  in  higher  plume  temperatures  and  a  delay  in  afterburning  shutdown.  This  may  be 
seen  in  Figs.  4  and  5  which  are  plots  of  C02  mole  fraction  and  temperature,  respectively,  at  x/xreji  =  3 
for  each  altitude.  The  transverse  coordinate  in  these  figures  has  been  scaled  by  5  which  is  the 
location  at  which  the  N2  mole  fraction  is  99%  of  its  ffeestream  value  for  the  laminar  rate  cases. 
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From  Fig.  4,  the  temperature  pdf  model  is  seen  to  enhance  the  formation  of  CO2  for  all  altitudes 
except  25  km.  At  this  altitude  the  laminar  rate  model  produces  more  CO 2  than  the  pdf  model. 
However,  the  H2O  formation  (not  shown)  at  25  km  is  much  greater  for  the  pdf  model,  offsetting  the 
lower  CO 2  formation  so  that  the  pdf  model  produces  higher  peak  temperatures  for  all  altitudes  (Fig. 
5).  The  enhancement  of  the  H2O  chain  within  the  mechanism  results  in  a  large  enhancement  of  OH 
for  the  pdf  cases  as  well  (Fig.  6).  For  the  laminar  rate  cases,  Fig.  4  shows  a  large  drop  in  peak  CO 2 
at  this  streamwise  location  as  afterburning  shuts  down  and  the  plume  ignition  point  moves  further 
downstream.  The  pdf  model  shows  a  similar  trend  but  to  a  lesser  extent,  indicating  the  plume  is 
shutting  down  slower.  Also,  similar  to  Fig.  3,  Fig.  5  shows  the  impact  of  the  pdf  model  on  the 
temperature  field  to  increase  with  increasing  altitude. 

The  large  impact  of  the  assumed  pdf  method  on  plume  combustion  processes  results  from  the 
influence  of  temperature  fluctuations  on  the  evaluation  of  the  mean  chemical  reaction  rate.  Fig.  7 
presents  a  plot  of  the  temperature  fluctuation  intensity  as  a  function  of  altitude  at  x/xreji  —  3.  The 
fluctuation  levels  are  seen  to  be  quite  large  at  this  streamwise  station  and  increase  in  magnitude  with 
increasing  altitude.  These  large  values  are  a  result  of  the  excitation  of  the  plume  shear  layer  at  the 
point  where  the  reflected  barrel  shock  intersects  the  layer.  The  increasing  magnitude  of  the 
fluctuations  indicates  the  pdf  model  is  becoming  increasingly  important  with  increasing  altitude. 
Consequently,  neglecting  these  fluctuations  would  be  inappropriate,  especially  as  afterburning 
shutdown  is  approached.  The  primary  effect  of  the  temperature  pdf  model  is  on  the  ignition  process 
so  that  prior  to  shutdown  the  model  is  much  less  significant  as  seen  for  the  25  km  case. 

From  Fig.  7  the  model  predicts  the  temperature  intensity  to  increase  with  increasing  altitude. 
From  Equ.  (12),  it  is  unclear  what  the  source  of  this  trend  is.  Fig.  8  presents  a  plot  of  the  production, 
dissipation  and  dilatation  terms  of  Equ.  (12)  at  x/xrefi  =  3  for  each  altitude.  From  this  figure  the 
magnitude  of  each  term  is  found  to  be  decreasing  with  increasing  altitude.  The  dilatation  term  also 


16 


has  both  positive  and  negative  values  and  is  small  compared  to  the  production  and  dissipation  terms. 
Consequently,  the  modeled  compressibility  term  in  Equ.  (12)  did  not  explicitly  play  a  significant  role 
in  the  evolution  of  the  energy  variance.  However,  compressibility  effects  are  implicitly  influencing 
Equ.  (12)  through  k  and  s  that  evolve  from  the  compressibility  corrected  k-e  model  used  in  the 
simulations.  Fig.  9  presents  the  sum  of  the  production,  dissipation  and  dilatation  terms  in  Equ.  (12) 
at  this  same  streamwise  location.  From  this  figure,  the  total  source  term  contribution  to  Equ.  (12)  is 
found  to  be  rapidly  decreasing  with  increasing  altitude.  This  indicates  that  the  increasing 
temperature  intensity  fluctuations  seen  in  Fig.  7  are  a  result  of  upstream  influences  and  not  a  result  of 
enhanced  relative  production  within  the  shear  layer.  Examination  of  contour  plots  of  turbulent 
kinetic  energy  and  temperature  intensity  show  the  plume  shear  layer  to  be  greatly  influenced  by  the 
reflected  barrel  shock.  This  reflected  shock  (see  Fig.  2)  excites  the  shear  layer  causing  the  layer  to 
thicken  and  it  enhances  the  production  of  turbulence  and  temperature  fluctuations.  This  trend  is 
consistent  with  the  work  of  Norris  and  Edwards,20  who  applied  a  much  more  general  unsteady  large- 
eddy  simulation  technique  to  compute  high-speed  reacting  exhaust  flows.  This  shock-shear  layer 
interaction  determines  the  turbulence  and  temperature  fluctuation  levels  which  feed  into  the  evolving 
shear  layer.  This  interaction  is  the  apparent  source  of  the  trend  observed  in  Fig.  7. 

The  large  increase  in  plume  temperatures  for  the  temperature  pdf  simulations 
correspondingly  results  in  a  large  increase  in  plume  radiative  emissions.  Figs.  10  and  11  present 
comparisons  of  station  radiation  and  total  radiant  intensity,  respectively.  From  Fig.  10,  the  station 
radiation  for  the  pdf  model  is  seen  to  be  significantly  larger  than  for  the  laminar  rate  model.  Burning 
for  the  pdf  model  at  35  km  is  also  evident  while  the  laminar  model  is  approaching  extinction.  The 
higher  station  radiation  predictions  in  Fig.  10  produce  significant  differences  in  the  total  plume 
intensity  as  seen  in  Fig.  11.  As  noted  earlier,  the  pdf  model  has  a  greater  influence  as  altitude  is 
increased.  This  trend  is  also  very  apparent  in  Fig.  11.  This  figure  also  shows  the  pdf  model  to  delay 
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the  onset  of  afterburning  shutdown  and  cause  the  plume  to  bum  to  higher  altitudes  than  for  the 
laminar  rate  model.  The  pdf  model  also  significantly  changes  the  afterburning  shutdown  rate, 
slowing  it  over  what  is  seen  for  the  laminar  rate  model. 

Though  changing  the  magnitude  of  the  emissions,  the  pdf  model  did  not  change’ the  character 
of  the  shutdown  event  seen  in  Fig.  11.  Both  the  laminar  and  pdf  models  show  a  gradual  type 
shutdown  behavior  which  is  characteristic  of  a  Damkohler  number  effect  as  described  by  Calhoon.4 
As  demonstrated  by  Calhoon,4  a  strain  rate  induced  extinction  phenomenon  may  be  the  source  of 
rapid  afterburning  shutdown  behavior  which  has  been  observed  for  some  missile  systems.  A  strain 
rate  extinction  mechanism  is  a  consequence  of  turbulence-chemistry  interactions.  However,  an 
extinction  mechanism  is  not  modeled  within  the  present  assumed  pdf  formulation.  Therefore,  the 
assumed  pdf  model  cannot  capture  such  a  phenomenon. 

As  mentioned  in  Sec.  II.,  the  coefficients  C/  and  C2  on  the  production  terms  in  Equs.  (12) 
and  (13)  were  specified  as  C;  -0.2  and  C2  =  0.5.  The  value  C2  was  specified  by  optimizing  the 
coefficient  to  match  the  scalar  fluctuation  data  of  Lockwood  and  Moneib21  for  a  low  speed 
nonreacting  jet.  Applying  the  same  optimization  procedure  for  Equ.  (12)  resulted  in  a  value  of  0.343 
for  Cj.  Figs.  12  and  13  present  the  sensitivity  of  the  temperature  and  temperature  intensity  with 
respect  to  variations  in  C/  at  x/xrefl  =  3  for  the  30  km  case.  From  Fig.  13,  the  temperature  intensity  is 
seen  to  vary  by  —  50  %>  over  the  range  0.1  3  C\  31  0.343  while  the  peak  temperature  in  Fig.  12  varies 
by  only  ~  6  %.  For  simulations  with  C/  =  0.343,  the  temperature  intensities  were  high  enough  to 
force  the  plume  afterburning  to  be  at  a  near  equilibrium  condition  for  each  altitude.  Under  this 
condition,  the  normalized  C02  and  H20  mole  fractions  were  invariant  with  altitude  and  afterburning 
shutdown  was  not  initiated  within  the  25  to  35  km  window  of  the  simulations.  This  result  was 
unrealistic  for  the  generic  booster  under  consideration.  This  also  indicates  that  the  low  speed 
calibration  is  inappropriate  for  application  to  the  high  speed,  highly  compressible  flows  of  the 
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present  context.  The  present  value  of  Cy  =  0.2  was  selected  because  it  was  close  to  the  calibrated 
value  and  allowed  afterburning  shutdown  to  be  initiated  within  the  25  to  35  km  window,  as  expected. 
The  lower  value  of  Cy  =0.1  showed  the  same  trends  for  mean  flow  quantities  as  for  the  present 
value  of  Ci  =  0.2,  but  to  a  slightly  lesser  extent.  The  value  selected  here  for  C /  was  deemed 
appropriate  since  the  objective  of  this  study  was  to  make  an  assessment  of  the  potential  impact  of 
turbulence-chemistry  interactions  on  exhaust  plume  simulations,  and  not  to  make  quantitative 
predictions. 

This  adjustment  of  Cy  from  the  low  speed  calibrated  value  to  the  smaller  value  of  0.2  is  also 
consistent  with  the  recent  large-eddy  simulation  (LES)  study  of  Calhoon,  et  al.  In  that  study,  LES 
was  used  to  investigate  the  effect  of  compressibility  on  temperature  fluctuations  in  planar  shear 
layers.  It  was  found  that  compressibility  effects  in  the  high  speed  flow  regime  substantially  reduced 
temperature  fluctuations  compared  with  what  was  found  for  low  speed  flows.  The  reduction  in 
magnitude  of  these  fluctuations  was  also  similar  to  what  have  been  accomplished  in  Fig.  13  by 
changing  Cy.  This  finding  suggests  that  the  modeled  compressibility  term  in  Equ.  (12)  is  inaccurate 
due  to  the  fact  that  it  made  little  contribution  to  the  evolution  of  the  energy  variance  as  seen  in  Fig.  8. 
This  reduction  of  Cy  has  in  some  sense  crudely  represented  the  effects  of  compressibility  as  found  by 
Calhoon,  et  al 22  However,  to  resolve  this  issue  will  require  the  development  of  an  accurate 
compressibility  correction  for  Equ.  (12)  using  LES  or  experimental  data.  In  addition,  a  more 
advanced  formulation  may  be  required  that  includes  a  length  scale  equation  for  the  temperature 
fluctuations.23 

To  investigate  the  effect  of  species  fluctuations,  Fig.  14  presents  a  comparison  of  the  mean 
temperature  at  x/xrefi  =  3  calculated  using  the  temperature  pdf  only  model  and  the  full  temperature 
and  species  pdf  formulation  described  in  Sec.  H  From  the  figure  it  is  apparent  that  the  full 
temperature-species  pdf  formulation  produced  virtually  identical  results  to  the  temperature  pdf  only 
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model.  This  was  found  to  be  true  for  the  other  flow  variables  as  well.  The  ineffectiveness  of  the 
assumed  species  pdf  model  to  significantly  change  the  results  can  be  traced  to  the  prediction  of  the 
scalar  energy,  0,  from  Equ.  (13).  Fig.  15  presents  a  plot  of  Q  across  the  plume  shear  layer  at  x/xrep  = 
3  calculated  using  Equ.  (13)  both  with  and  without  the  chemical  source  term  at  30  km. '  Including  the 
chemical  source  term  causes  O  to  be  driven  to  small  values  effectively  turning  off  the  species  pdf 
model,  producing  the  same  result  as  for  the  temperature  pdf  only  model.  Eliminating  this  source 
term  from  Equ.  (13)  produces  large  values  of  scalar  energy  as  seen  in  the  figure.  This  chemical 
source  term  was  found  to  be  dissipative  and  to  destroy  the  scalar  energy  and  drive  it  to  small  values. 
Baurle,  et  al24,  observed  the  same  behavior  for  a  high  speed  H2-air  flame.  In  another  publication, 
Baurle,  et  al?5,  also  compared  predictions  of  this  source  term  using  the  assumed  pdf  model  and  a 
more  comprehensive  pdf  evolution  equation  model  for  the  same  jet  flame.  This  comparison  showed 
the  predictions  of  this  term  to  generally  be  of  opposite  sign  and  to  be  different  in  magnitude  by  a 
factor  of  five.  The  present  results  support  the  assertion  of  Baurle,  et  al 25 ,  that  the  present  form  of 
the  species  pdf,  Pq,  is  incapable  of  computing  the  higher  order  moments  required  to  evaluate  the  Q 
chemical  source  term  with  any  reasonable  degree  of  accuracy.  This  also  casts  doubt  on  the  assumed 
species  pdf  model’s  accuracy  for  computing  the  mean  species  reaction  rate  (Equation  (1))  which 
often  contains  higher  moments  as  well.  This  deficiency  of  the  model  is  likely  the  source  of  the 
ineffectiveness  of  the  assumed  species  pdf  formulation. 

As  a  final  topic,  to  demonstrate  the  adequacy  of  the  computational  grid  resolution,  Figs.  16 
and  17  present  plots  of  the  temperature  and  CO 2  mole  fraction  at  two  locations  in  the  plume  for  each 
altitude  condition.  Shown  in  these  figures  are  results  for  the  plume  solutions  on  the  original  (fine) 
grid  and  a  coarse  grid  with  half  as  many  points  in  each  coordinate  direction.  As  seen  in  Figs.  16  and 
17,  the  results  are  virtually  identical,  indicating  that  the  original  grids  adequately  resolved  the 
flowfields. 
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VI.  CONCLUSIONS 


A  computational  study  was  undertaken  to  assess  the  impact  of  turbulence-chemistry 
interactions  on  the  afterburning  and  afterburning  shutdown  characteristics  of  a  generic  amine 
booster.  Turbulence-chemistry  interactions  were  accounted  for  using  an  assumed  probability  density 
function  based  method.  Analysis  of  the  simulation  results  including  the  model  led  to  the  following 
conclusions: 

1)  Turbulence-chemistry  interactions  modeled  using  the  assumed  pdf  method  were  found  to 
be  a  first  order  effect  and  have  a  large  impact  on  plume  signatures  as  afterburning  shutdown  was 
approached.  The  pdf  model  enhanced  plume  afterburning  and  delayed  the  onset  of  shutdown  causing 
the  plume  to  bum  to  higher  altitudes  than  when  the  model  was  not  included.  However,  the  pdf 
model  did  not  change  the  character  of  the  shutdown  event.  The  model  produced  a  gradual  type  of 
shutdown  behavior  characteristic  of  a  Damkohler  number  shutdown  mechanism  (Calhoon,  2000).4 
The  model  does  not  account  for  the  rapid  shutdown  phenomenon  identified  by  Calhoon.4 

2)  The  assumed  pdf  model  was  found  to  become  increasingly  important  as  altitude  increased. 
This  resulted  from  the  increasing  trend  of  temperature  fluctuations  within  the  shear  layer.  This  trend 
was  found  to  be  a  result  of  shock-shear  layer  interactions  within  the  plume. 

3)  The  calibration  constant  on  the  production  term  in  the  energy  variance  equation  obtained 
from  optimization  using  low  speed  jet  data  was  found  to  be  inapplicable  to  high  speed  plume  flows. 
This  low  speed  calibration  produced  excessively  high  temperature  fluctuations  that  resulted  in 
unrealistic  results.  Lower  values  of  the  coefficient  produced  realistic  results  and  demonstrated  the 
importance  of  turbulence-chemistry  interactions  within  plume  flows.  A  recent  LES  study22  suggests 
that  the  source  of  the  difficulty  is  the  lack  of  a  credible  compressibility  correction  within  the  energy 
variance  equations.  Consequently,  a  reduction  in  the  production  calibration  constant  was  required  to 
bring  the  modeled  scalar  fluctuation  statistics  in  line  with  LES  data.  This  issue  highlights  the  need 
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for  a  rigorous  validation  of  the  energy  variance  equation  for  high  speed  flows  before  quantitative 
predictions  can  be  made. 

4)  The  species  fluctuation  aspect  of  the  assumed  pdf  formulation  was  found  to  be  ineffective 
in  the  present  plume  flows.  Apparent  inaccuracies  in  the  evaluation  of  the  chemical  Source  term  in 
the  scalar  energy  equation  resulted  in  the  model  turning  itself  off.  This  same  trend  has  been 
observed  by  other  investigators  and  suggests  that  a  reformulation  of  the  species  fluctuation  aspect  of 
the  model  is  warranted. 
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Fig.  1.  Illustration  of  the  behavior  of  plume  total  radiant  intensity  as  a  function  of  altitude,  (a) 
afterburning  regime,  (b)  afterburning  shutdown  regime,  and  (c)  post-afterburning  shutdown  regime. 
Fig.  2.  Contours  of  (a)  temperature  and  (b)  CO 2  mole  fraction  at  30  km  for  both  the  laminar  rate 
(lam)  and  temperature  pdf  (t  pdf)  models.  The  nozzle  exit  plane  is  at  the  top  and  downstream  is  to 
the  bottom.  The  transverse  direction  has  been  scaled  by  a  factor  of  5  for  clarity. 

Fig.  3.  Comparison  of  laminar  and  temperature  pdf  rate  model  predictions  of  centerline  mean 
temperature  at  25, 30  and  35  km  (x  =  0  is  at  the  engine  nozzle  exit  plane). 

Fig.  4.  Prediction  of  transverse  mean  C02  mole  fraction  at  x/xre/i  =  3  as  a  function  of  altitude  for  the 
laminar  rate  and  temperature  pdf  models. 

Fig.  5.  Prediction  of  transverse  mean  temperature  at  x/xreji  =  3  as  a  function  of  altitude  for  the 
laminar  rate  and  temperature  pdf  models. 

Fig.  6.  Prediction  of  transverse  mean  OH  mole  fraction  at  x/xrefl  =  3  as  a  function  of  altitude  for  the 
laminar  rate  and  temperature  pdf  models. 

Fig.  7.  Temperature  intensity  predictions  from  the  assumed  pdf  model  as  a  function  of  altitude  at 
x/x^eji  3. 

Fig.  8.  Transverse  variation  of  the  production,  dissipation  and  dilatation  terms  in  Equ.  (12)  at  x/xrefl 
=  3  as  a  function  of  altitude. 

Fig.  9.  Transverse  variation  of  the  sum  total  of  the  production,  dissipation  and  dilatation  terms  in 
Equ.  (12)  at  x/xreji  =  3  as  a  function  of  altitude. 

Fig.  10.  Source  station  radiation  variation  as  a  function  of  altitude  at  a  90  degree  aspect  angle. 

Fig.  11.  Comparison  of  laminar  rate  and  temperature  pdf  model  predictions  of  source  total  radiant 
intensity  as  a  function  of  altitude  at  a  90  degree  aspect  angle. 

Fig.  12.  Sensitivity  of  mean  temperature  to  the  calibration  constant  C}  in  Equ.  (12)  at  x/xrefl  =  3  for 
30  km. 
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Fig.  13.  Sensitivity  of  temperature  intensity  to  the  calibration  constant  C;  in  Equ.  (12)  at  x/xrejj  =  3 
for  30  bn. 

Fig.  14.  Prediction  of  transverse  mean  temperature  at  x/xreji  =  3  as  a  function  of  altitude  using  the 
temperature  pdf  only  (t  pdf)  and  the  temperature-species  pdf  (ty  pdf)  formulations. 

Fig.  15.  Prediction  of  scalar  energy  from  Equ.  (13)  both  with  and  without  the  chemical  source  term 
included  at  x/xrejj  =  3  and  30  bn. 

Fig.  16.  Plume  temperature  profiles  for  the  original  (fine)  and  coarse  grid  solutions  for  the  laminar 
£ 

rate  case  at  x  =  x/xre/i  =  3  and  9. 

Fig.  17.  Plume  CO 2  mole  fraction  profiles  for  the  original  (fine)  and  coarse  grid  solutions  for  the 
laminar  rate  case  at  x*  =  x/xreji  =  3  and  9. 
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